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Abstract 

We associate with k hermitian TV x TV matrices a probability measure on M. k . It is 
supported on the joint numerical range of the fc-tuple of matrices. We call this measure 
the joint numerical shadow of these matrices. Let k = 2. A pair of hermitian TV x TV 
matrices defines a complex TV x TV matrix. The joint numerical range and the joint 
numerical shadow of the pair of hermitian matrices coincide with the numerical range 
and the numerical shadow, respectively, of this complex matrix. We study relationships 
between the dynamics of quantum maps on the set of quantum states, on one hand, and 
the numerical ranges, on the other hand. In particular, we show that under the identity 
resolution assumption on Kraus operators defining the quantum map, the dynamics 
shrinks numerical ranges. 

Keywords: joint numerical range, positive matrices, quantum maps, quantum states 
2000 MSC: 47A12, 81P16, 51M15 



1. Introduction 

We denote by Hn the space <C N endowed with the standard scalar product (<p\ij)). 
Thus, "H^r is the Hilbert space of TV dimensions. Let A be an operator on Hn- Its 
numerical range (also called the field of values) W(A) C C is the set of numbers z = 
(ip\A\tp), where (iplip) = 1 QQ- A crucial fact, conjectured by O. Toeplitz in 1918 Q 
and proved by F. Hausdorff in 1919 0] is that W(A) is convex. See 0] for a modern 
exposition and historical remarks. If A is a normal operator, then W{A) is the convex 
hull of the spectrum of A [l[, thus a convex polygon. For generic A the boundary dW(A) 
is smooth 0. For TV = 2, dW(A) is a (possibly degenerate) ellipse 0,0]. See Q for a 
classification of numerical ranges when TV = 3. 

In recent years, numerical range and its generalizations found novel applications in 
quantum mechanics, in particular, in the theory of quantum information Hol . 11 1. 



We will freely use the physics terminology in what follows. Pure quantum states are 
the hermitian projections of Hat onto one-dimensional subspaces C| - 0). The set J7jv of 
pure quantum states is naturally isomorphic to the complex projective space CP W_1 . A 
density matrix is a positive operator p > on H.^ satisfying Trp = 1. A pure quantum 
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state is a density matrix, and any density matrix is a convex combination of pure quantum 
states. In the physics terminology, density matrices are the quantum states. The set Qn 
of quantum states is convex, and f2jv C Qn is the set of its extremal points. The elements 
of Qn \&n are mixed quantum states. 

Applications of the numerical range to quantum information are based on the inter- 
pretation of as a plane projection of £In [Hj]. For instance, in the one qubit case, 
A = 2, we have Q 2 = CP 1 = S 2 (the Block sphere), and Q 2 = B 3 C R 3 (the Block ball). 
A projection of the two-sphere is a possibly degenerate ellipse. This simple observation 
is the base for the well known facts about numerical ranges of 2 x 2 matrices Q • 

Let A be an operator on Hn- Its numerical shadow Pa(z) is a probability distribution 
on C supported on W(A) [H 03, El- Let be distributed on the unit sphere in 
Hn according to the unique rotationally invariant probability measure (i.e., the Haar 
measure). Then Pa{z) is the probability density of z — (ij)\A\ip). Let A = 2. Then 
Pa{z) is the density of the plane shadow of the sphere under a light beam. This explains 
the term "numerical shadow" (lBj . 

Let Ai,.. ■ A m be hermitian operators on Hn- Their joint numerical range (JNR) 



W(Ai, . . . , A m ) c R m [ljfl is defined as follows: 



W(A 1 ,...,A m ) = {((^|Ai|^), . . . , {ip\A m \ijj) ) : (ip\ip) = 1}. (1) 

Let m = 2. Then W(Ai,A 2 ) = W(A\ + iA 2 ). Thus, JNR generalizes the notion of the 
numerical range of a complex operator. If m > 2, then W(A\, . . . ,A m ) C R m is not 
necessarily convex For instance, JNR of the triple of Pauli matrices is the Bloch 
sphere. We will elaborate on this observation in section [2] 

In this work we establish and explore a relationship between the notion of the joint 
numerical range and the set of quantum states. In Theorem [T]we show that the JNR of an 
?7i-tuple of hermitian operators is a suitable linear projection of the set of pure quantum 
states. Using this observation, we associate with the joint numerical range a probability 
measure on R m , the joint numerical shadow of an m-tuple of hermitian operators. See 
section This extends the concept of numerical shadow to arbitrary dimensions (TBI . 
13 . Il4 . 1 X 61 ] . Let $ denote the discrete dynamics on the set of quantum states defined by a 
set of k Kraus operators satisfying the identity resolution equation (fTTj) . Wc analyze the 
effect of the dynamics $ on the numerical range. In Corollary [TJwe obtain an inclusion 
of the numerical range of $(A) into the numerical range of A. 



2. Joint numerical ranges 

Let £jv (resp. Mn, Vn, Qn) denote the space of all (resp. hermitian, positive 
definite, positive definite with trace 1) linear operators on T-Ln- Let Q,n C Qn be the set 
of rank one projections. Then Cn endowed with the scalar product 

(A, B) = tr(AB*) (2) 

is a Hilbcrt space of A 2 dimensions, and A4n C Cn is a real subspace of A 2 dimensions 
on which the scalar product @ induces the structure of a euclidean space. The set 
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Vn C A4n is a closed convex cone. Its interior consists of strictly positive operators, 
p > 0, and its apex is the zero operator. The set Qn is the intersection of Vn and the 
hyperplane {Trp = 1}. It is a bounded convex region (i.e., has nonempty interior) in the 
(A^ 2 — l)-dimensional affine space, and f2jy C Qn is the set of its extremal points. For 
a unit vector ip £ TLn we set p^ = \ip)(ip\ £ fijv- The manifold <CP N ~ 1 is the quotient 
of the unit sphere S(Hn) by the natural linear action of the unit circle S 1 C C. The 
map iJj i — y Pip yields an isomorphism of CP W_1 and In what follows we will identify 
CP N ~ X and £In via this isomorphism. 

Let Ai, . . . ,A m £ A4n be arbitrary. The mapping from CP N_1 to R m given by 

j™ 4, 4, : M MAity), . . . , (V|A m |^)). (3) 

is the jomi numerical range map. The range W(.Ai, . . . , A m ) of this map is the joint 
numerical range (JNR) of operators Ai, . . . , A m . By the isomorphism t CP N ~ 1 ~ f2jvj we 
have j nr (Ai, ...,A m ) '■ 

We will recall a few basic notions in affine geometry. By a vector space we will mean 
a finite dimensional real vector space. Let V be a vector space. A set H C V is an affine 
subspace if there is a linear subspace Hq C V and a vector ho £ H such that 77 = Ho + ho. 
Let G C J7, H C V be affine subspaces. Let GoC^ffoC^ be the corresponding linear 
subspaccs. A map Af : G — > H is an affine isomorphism if there is a linear isomorphism 
A fo ■ Gq —> Ho and vectors ho £ H, go £ G so that 

Af(g) = Afo(g-g ) + h . 



Definition 1. Let U, V be vector spaces, and let X C U,Y £\V be arbitrary sets. They 
are affincly isomorphic if there exist affine subspaces H C U,G C V such that X C 
G,Y C H , and an affine isomorphism A : G — > H such that A{X) = Y . The induced 
map A : X — > Y is an affine isomorphism of X onto Y. 

Note that linear isomorphism of sets are the special cases in the above setting when the 
subspaces and the maps in question are, actually, linear. We will not distinguish between 
the linear and affine situations in what follows. 

We will now expose a general topic in linear algebra. Let U, V be vector spaces. We 
assume that U is a euclidean space with the scalar product (•,•). Let m > 1. With 
any nonzero vectors u±, . . . , u m £ U and vx, . . . , v m £ V wc associate a linear operator 
L : U -> V by 

m 

L ( u ) = ^2(u,Ui)vi. (4) 

i=i 

Let A C U (resp. B C V) be the subspace spanned by u%, . . . ,u m (resp. v±, . . . ,v m ). 
We will need a simple lemma about the operator in equation (U). The following claim 
should be in the literature. We do not give a reference, but the proof is straightforward, 
and wc leave it to the reader. 
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Lemma 1. 1. Let dim .4 = p < m. Assume without loss of generality that the vectors 
Uy, . . . , Up span A. Then there are vectors Vy, . . . , v p S B such that the linear operator in 
equation ^ satisfies 

v 

L(u) = y](u,Ui)vj. (5) 

i=l 

2. Let dimi3 = q < m. Assume without loss of generality that the vectors vy, . . . ,v q span 
B. Then there are vectors Uy,...,u q € A such that the linear operator in equation (j4]) 
satisfies 

<? 

L(u) = ^(u, Ui)vi. (6) 

i=l 

3. If p = m (resp. q = m) then Vi = Vi (resp. ui = Ui) for 1 < i < m. 

The proposition below is immediate from Lemma [TJ 

Proposition 1. Let the setting be as in Lemma [TJ and let Pr_\ : U — > A be the or- 
thogonal projection. Then i) There is a unique linear operator M : A — > B such that 
L = M o Prj^; ii) If dirndl = m (resp. dimS = m) then M is surjective (resp. infec- 
tive); Hi) If dirndl = dim£? = m then M is an isomorphism. 

Let G, H C U be affine subspaces in a vector space. We say that they are parallel if 
H = G + uq for some uq £ U. 

Corollary 1. Let U,V be vector spaces, and let L : U — > V be as in equation ([¥]). 
Let A C U, B C V be the subspaces spanned by the vectors u\, . . . , u m and v\, . . . , v m 
respectively. Let M : A ^ B be the operator from Proposition [7J 

Let G C U be an affine subspace containing a subspace parallel to A. Let V C G be an 
arbitrary set. If M is injective, then Pr^(T) and HT) are affinely isomorphic. Proof. 
Let uq € U be such that A = G + uq. Let Bo (Z B be the range of M . Thus, M : A — > Bo 
is a linear isomorphism. By Proposition [7J we have 

L(T + u Q ) = M Pr A {T + u ) = M(T + u ) C B Q . 

The mapping u <— > L(u + Uq) induces an affine isomorphism of G and Bq. ■ 

We will now apply the above material to the joint numerical range. 

Theorem 1. Let A\, . . . , A m G Mn be traceless, linearly independent hermitian oper- 
ators. Let A C Mn be the subspace spanned by them. Then i) The joint numerical 
range of Ay, . . . ,A m is affinely isomorphic to Pr^(flN); H) The convex hull of the joint 
numerical range of Ay, ... , A m is affinely isomorphic to Pr^Q^)- 

Proof. Let U = Mn, V = W" 1 . Let Vy, . . . ,v m be the standard basis in W l . Then the 
map jnrA 1 ,...,A m has the form equation (0]) with Ui = Ai. By Proposition [TJ the map 
M : A M m 'is a linear isomorphism. Set G = {X e M N ■ tr(X) = 1}. Then G C M N 
is an affine hyperplane containing the affine subspace A+ In- Claim i) now follows from 
Corollary [T] 

The convex hulls of afinely isomorphic sets are afinely isomorphic. Hence, claim i) 
implies claim ii). ■ 
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Corollary 2. Let the setting be as in Theorem Q with m = N 2 -l. Then W(A ll . . . , A m ) 
is affinely isomorphic to Q,n and conv[W(Ai, . . . , A m )] is affinely isomorphic to Qn ■ 
Proof. This is a special case of Theorem [7J In this case the affine hyperplane G is 
parallel to A. Hence Pr^ induces an affine isomorphism of fijy an d Pta{^n)- ■ 

Example 1. Let N = 2. The Pauli matrices 
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' -i ' 
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0"! = 


1 





i ^2 = 
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, ^3 = 


-1 



(7) 



form an orthogonal basis in the space of traceless Hcrmitian 2x2 matrices. By Corol- 
lary [2^1 W{a\, 02, o~s) C R 3 is a sphere. Let p be a 2 x 2 density matrix. Its Bloch vector 
representation is given by 

1 3 

j'=i 

We have Tj = tr(ojp); the Bloch vector (n, T2, 73) = r G R 3 determines p. The positivity 
condition p > 0, implies 

\Jt?+t1 + t$ = \\t\\<1/2. 

The equality ||r|| = 1/2 holds iff p = p^,, a pure state. Thus, the map j nr (o-i,cr 2 ,cr 3 ) yields 
an isomorphism of Q2 (resp. ^2) and the Bloch ball (resp. Bloch sphere). 

We will now discuss a generalization of Example 1 to N > 2. Note that d = N 2 — 1 
is the dimension of the space of traceles N x TV hermitian matrices endowed with the 
scalar product equation ((3]). Let Ai, . . . , Ad be a basis in this space. It is convenient to 
use an orthogonal, but not necessarily an orthonormal basis. For example, for N = 3 it 
is standard to use the Gell-Mann matrices (Ai, . . . , As). They satisfy the orthogonality 
relations tr(AjAk) = 2<5jk 0. Let p be a N x N state. The counterpart of equation © 
is 

1 d 

The generalized Bloch vector r = r(p) has d components Tj = Tr(pAj)/Tr(A^). Assume 
now that the basis Ai,...,A<j is orthonormal. Then the generalized Bloch vector 
of a pure state p^, = \ip){tp\ satisfies tj = (^|Aj|V')' Hence the joint numerical range 
W(Xi, . . . Ad) contains the Bloch vector r^. 

By Corollary [5J the convex hull conv[W(Ai, . . . Ad)] is affinely isomorphic to the set 
Qn of quantum states. By Theorem [TJ if Ai, . . . , A m are linearly independent, then the 
joint numerical range W(Ai, . . . , A m ) is affinely isomorphic to a projection of Qn to R m . 
This is a generalization of a result in [l2[ . 

3. Quantum maps 

We recall that the set Cn of operators on Hn is a Hilbcrt space with the scalar 
product (A, B) = tr(AB*). Let CCn be the space of linear operators on Cn- For 



1 See also [HEl. 
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$ € CCm we denote by <&* S ££at its adjoint. For X, Y e £yv we define <&x|y € CCm 
by 

$x|y(A) = XAY*. (10) 

The following lemma will be used in the proof of Proposition [3] The proof is straight- 
forward and we leave it to the reader. 

Lemma 2. 1. The operators §x\y span the vector space CCm- 

2. Let X 1 ,Y 1 ,X 2 ,Y 2 e C N . Then 

®x 2 \y 2 ®x 1 \y 1 = ®X 2 X 1 \Y 2 Y 1 - 

3. We have 

We will denote by 1 <G Cm the identity operator. 

Definition 2. Let k > 1. We say that the k-tuple of operators X\, . . . ,X^ € Cm is an 
identity resolution if 

k 

J2XiX* = l, (11) 

i=l 

The dual property 

k 

J2x*Xi = l, (12) 

i=l 

holds iff XI , . . . , X£ is an identity resolution. Let Xi,..., Xk G Cm be arbitrary. Set 

h 

®{x u ...,x k )=^$x i \x i - (13) 

i=l 

The following lemma summarizes the basic properties of operators 3>(Xi,....x fc ) ■ 

Lemma 3. Let X\,...,Xk S Cm be arbitrary, and let <f> = §(Xi,...,Xk) € CCm- Then 
the following holds. 

1. The operator $ : Cm Cm preserves the sets M.m ond Vm ■ 

2. The operators X\, . . . , Xk satisfy equation (fTTj) iff $ : Cm — > Cm preserves 1. 

3. The operators X\, . . . ,Xk satisfy equation (|12p iff $ : Cm — > £./v preserves the trace 
iff <!>* preserves 1. 

^. If X\, . . . ,Xk satisfy equation (fl"2"l) . i/ien $ : £at — > £jv preserves Qm ■ 
Proof. It suffices to prove 1 ) for operators $ = <&x\x ■ The former property is immediate 
from (XAX*)* = XA*X*, and the latter from (XAX*v,v) = (AX*v,X*v). Claim 2) 
is immediate from the definition o/$. We have 

TV#(A) = ($(A),1) = (A,$*(1)). 

Claim 3) now follows from Lemma [H Claim 4) is immediate from 1) and 3). ■ 
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For tjj G Hn we set p^ = \tp)(ip\. Then p^ G Vn, and for any a G C we have 

pai> = M%- (14) 
In particular, p^ = iff tp = and p,/, G Qat iff ||V'|| = 1- 

Proposition 2. Lei Xi,...,Xk G £jv, arcd /e£ $ = Let -0 G "Htv &e any 

vector. For 1 < i < k such that Xiip ^ 0, set ipi — Xiip/\\Xiip\\. 

Then 

*0V) = XI||^||V- (15) 
Proof. Let 1,7 G £_/v and u G Hat- TTien 

(*x\yM)v = X{^,Y*v)i> = {Y^,v)Xi>. 

In particular, we have 

®x\x{Pi>) = Px^- (16) 
The claim now follows from equations (|13|) and (1141) . ■ 



Corollary 3. LetXi, . . . ,Xk satisfy equation (JT2J) andsetQ = $(Xi,...,x k ) ■ Lei V G Hn 
be a unit vector. Let 1 < m < k be the number of unit vectors ipi from Proposition [Jl 
Then 

m 

= 'Yl/PiP^n 

i=l 

where p\ , . . . , p TO > and = !• 

Proof. Set pi = \ \Xirt>\\ 2 . The claim follows from equation (JTSJ) and t/ie identity 



Corollary 1. Let Jfj, . . . G £jv ; set $ = Suppose that X\, . . . ,Xf. 

satisfy equation (|11[) . TTien 
L For any A G £tv we have 

W(Q(A)) C (17) 
For any Ai , . . . , A m G M. n we have 

W($(Ai), . . . , $(A m )) c conv(W(Ai, . . . , A m )); (18) 

3. If W(Ai, . . . , A m ) is convex, then we have 

W(*CAi), . . . , $(A m )) c W(i4!, . . . , A m ). (19) 
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Proof. Let A <G £n be arbitrary, and let i/j £'Hn be a unit vector. Then 

Tr($(A) PV) ) = Tr (^^AX*)^ = Tr UQ^Xf p^X{)\ . 

Note that the operators Aj* , . . . , X£ satisfy equation (|12p . By Proposition and Corol- 
lary [3 there are unit vectors ipj € Hn and probabilities pj such that 

Tr(4>(A)^) =^ Pj Tr(A^) = ^ Pj (^|A|^>. (20) 
j j 

Claim 2 now follows from equation ^ . Claim 3 is a special case of Claim 2. Claim 3 
and the Toeplitz-Hausdorjf theorem yield Claim 1. ■ 



Let $ = $ 



(*i 



X k )- In the physics literature the transformations of Vn defined by 



(21) 



4=1 



are called quantum maps. Operators Xi £ £n are the measurement operators or Kraus 
operators. Quantum maps correspond to generalized quantum measurements with k 
possible outcomes [l8j |. By Lemma [3j the map $ : p i— > p± is unital (resp. irace 
preserving) iff equation (|11[) (resp. equation f|12|) ^ is satisfied. The quantum map 



<[>* = $ 



x*) is dual to <&. The duality of quantum maps corresponds to the duality 



between the Schrodinger and the Heisenberg representations in quantum mechanics. In 
the former, the quantum states p G Vn evolve via p\ = <&(p), while the observables 
A E A4n do not. In the latter, the quantum states do not change and the observables 
evolve by A i-> A± = The identity 

Tr($(p)A) = Tr (s^X^a) = Tr {^ P X*AxA = Tv{p^{A)) 

illustrates this duality. Note that <J> is trace preserving (resp. unital) iff $* is unital 
(resp. trace preserving). Corrolary [1] describes a relationship between quantum maps 
and numerical ranges. We will now present a few examples illustrating this relationship. 



Example 2. Let N — 2. The decaying channel is the discrete dynamics $ 



corresponding to the Kraus operators X\ 



1 







and X? 



VT^P 

p G (0, 1) is a free parameter. Set = $ J_1 (A (1) ), where A 1 - 1 



Vp 


1 




, where 
Note that 



tr(A' 1 ^) = 0. The set W(A^) is the disc of radius 1/2 centered at the origin, i.e., 
at the baricenter point tr(AW)/2. We have A® = (1 - p)^'" 1 )/ 2 ^ 1 ), and W(A^) 
: 1 - p)^" 1 )/ 2 ^^ 1 )) is the disc of radius ±(1 - p)^- 1 )/ 2 with the center at 0. For 

and A® = &(AW) = 



instance, A^ = 









1-p 





The 
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-0.25 
-0.5O 



50-0.25 0.00 0.25 0.50 
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0.50 
0.25 
O 0.00 
-0.25 
-0.50 



.50-0.25 0.00 0.25 0.50 



Figure 1: Numerical ranges W(AW), WiA^), W(A^) from Example 2 for p = 0.5 and 
W{BW), VK(B (2) ), W(B^) from Example 3 for p = 0.25. The stars mark the baricenter points in 
these examples. They are preserved by the dynamics. 



limit of W(AW), as j -> oo, is {0}. Fig. la shows W(A^), W{A&), W(A^) for 
p = 0.5. 

Example 3. Let again N = 2. The phase-flip channel is the discrete dynamics \E' = 
^x lt x 2 corresponding to the Kraus operators X\ = y/1 — p 1 2 and X 2 — y/pcri. Here 
again, p G (0, 1) is a free parameter. We set = A^ and fiW -1 ) = Then 

1 - 2p(l - p) " 
2p(l-p) 



B (2) 



1-p 
p 



and £ (3) 



. The numerical ranges of 



tors are X 1 = y/l—p — ql 3 , X 2 = sfp 



all are ellipses. Fig. lb shows W(B^),W(B^),W(B^) for p = 0.25. 

Example 4. This example is a generalization of Example 3 to N = 3. The double flip 
channel acting on a qutrit is the the discrete dynamics S = §x 1 ,x 2 ,x 3 ', the Kraus opera- 

oioi r i o o " 

10 and X 3 = y/q 1 .The pa- 

1 J L 1 

rameters p, q satisfy < p, q, p+q < 1. The operators X 2 , X% correspond to bit flips with 

[01 

probabilities p and q. The operator S is trace preserving. Let C^ 1 ' =010. This 

2i 

operator was studied in flij : its numerical range W{C^) is formed by the convex hull 
of an ellipse and a point outside it. The baricenter of W (C^) is tr(CW)/3 = (l + 2i)/3. 

p 1-p-q q 

Set = E^-^C^). For instance, = p 1 -p- q(l - 2i) 

2i + q{l - 2i) _ 

In Figure [Hp = 0.5 and q = 0.4, hence 1 — p — q = 0.1 and all three Kraus operators 
contribute to the dynamics. 
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-0.5 0.0 0.5 1.0 1.5 

Figure 2: This figure illustrates Example 4 and the inclusion relation equation 117t . It shows the 
numerical ranges W{C^) D W(C^~>) D W{C { '^) for p = 0.5, q = 0.4. The limit of W{C^), as j -> oo, 
is the baricenter point (1 + 2£)/3 marked by the star in the figure. 



4. Joint numerical shadows 

The set S(H N ) = {i> £ Wjv : ||V>|| = 1} is the unit sphere S 2N ~ l . We denote by 
the normalized Haar measure on S(Hn), as well as its push-forward to fl^ — CP W_1 
Let va be the push-forward of fx to C under the numerical range map of an operator 



A £ Cn- Then va is the numerical shadow of A [14l. 1 1 3l| . Let dz denote the density of 
Lebesgue measure on C. Then 4va(z) = PA(z)dz where the probability distribution is 
given by 

Pa(z)= [ (22) 

J O 

Let now Ai, . . . ,A m £ Mn, let jnr/^ A \ : f2jv — > R m be the joint numerical range 
map, and let W^-Ai, . . . , A m ) be the JNR. The push-forward v Alt ..., Am of [i to W % under 
the map jrrr^ A \ is the joint numerical shadow of A\, . . . ,A m . Thus, v Al ,,..,A m is 
supported on W(.Ai, . . . , A m ). The corresponding probability distribution satisfies 

m 

P Au ...,A m {x 1 ,...,x m )= duty) l[6(x j -( 1 p\A j \ip)Y (23) 

Numerical shadow is the special case of the joint numerical shadow corresponding to 
m = 2. See the examples below for illustration. 

If v is a measure on K m and a £ R, we denote by v a the push-forward of v under 
the self-map v i-> av of R m . If v\ , ^2 are measures on R m , we denote by v\ * v^ their 
convolution. The following proposition exposes a few basic properties of joint numerical 
shadows. We leave the proof to the reader. 



2 It is known as the Fubini— Study measure in the physics literature. 
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Proposition 3. 1. Let A\, . . . ,A m <E Mn- Let U G £n be a unitary operator. Set 
Bi = UAiU*, l<i<m. Then 

VB!,...,B m — VA ll ...,A m - 

2. Let A\, . . . , A m € Mn and B\, . . . , B m € Mn be arbitrary. Let a, b G K. Then the 
joint numerical shadows satisfy 

VaA 1 +bB 1 ,...,aA m +bB m = ^A 1 ,...,A m * v B u ...,B m - 



Example 5. We review Example 1. The set W(ai, 0-2,0-3) C M 3 is the sphere of radius 
1/2. The joint numerical shadow i> ai . CT2 ,o- 3 is the normalized Haar measure. 

Example 6. We use the isomorphism H4 = H2 ® H-i to define the extensions of Pauli 
matrices Aj = o-j ® I2, j = 1, 2, 3. We use the fact that the Haar measure on the space of 
pure states in T-La ® Wb induces by partial trace, lu = Trs|^i}(^|, the Lebesgue measure 
on the space of mixed states in [201 ] . Using the equality between the expected values 
of an operator on and the extended operator on H4, we obtain that VAx.A2.A3 is the 
normalized Lebesgue measure on the Bloch ball. Let Bj 
operator 



I2 ® o-j, j = 1, 2, 3. The swap 



(24) 



is a unitary operator on T-L^ satisfying Bj = SAjS*. By Proposition [5] and the above 
discussion, vbi,b 2 ,b 3 is the normalized Lebesgue measure on the Bloch ball. 
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